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into account. In this model, the general form of the susceptibility and group velocity of 
the probe laser we obtained are operators. Their expectation value and fluctuation can 
be obtained on the Fock space. Furthermore the uncertainty of the group velocity under 
very weak intensity of the controlling laser and the uncertainty relation between the phase 
operator of coupling laser and the group velocity are approximately given. Considering 
the decay rates of various levels, we may analyze the probe laser near resonance in detail 
| and calculate the fluctuation in both absorption and dispersion. We also discuss how the 

fully quantized model reduces to a semiclassical model when the mean photon numbers 
of the coupling laser is getting large. 
PACS numbers: 42.50.Gy, 42.65.-k, 42.65.An 
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Controlling the phase coherence in ensembles of multilevel atoms has led to the observation 
of many striking phenomena in the propagation of near-resonant light. The notable examples 
include ultraslow light pulse propagation, light storage, lasing without inversion and EIT[1, 
2, 3, 4, 5, 6]. Among these striking phenomena, ultraslow light speed and superluminality 
look more attracting. The drastic reductions in the group velocity of pulses were discussed 
in Ref.[5] and recent experiments have taken the reduction of the speed of light to extreme 
limits (shown in Refs.[6, 7] and even to zero[8, 9]. On the other way, the superluminal 
velocity was also observed in the abnormal-dispersion media[10]. In the most previous work 
the phenomenon of EIT and the accompanying enhancement of the index of refraction and 
susceptibility are treated using semiclassical theory in which both the coupling and probe 
lasers were treated as classical. In such a treatment, the occurrence of EIT requires the 
coupling laser to be much stronger than the probe laser. In Ref.[8], Fleischhauer and Lukin 
treated the probe laser as quantized, and showed that the quantum description of laser is 
more fundamental than the classical one, having advantages in uncovering new effects of EIT. 
In Ref.fll] Kuang et al. treated both the coupling and probe lasers as quantized and they 
pointed out that in general the group velocity depends on the intensity of the coupling as well 
as probe laser. But they ignored decay rates of various levels, so their treatment is essentially 
a time-independent approach. Obviously, if the decay rates are incorporated, it needs to solve 
the evolution equation of the density matrix, which is difficult in a fully quantized treatment. 

Based on the results of the previous works, the strength of the coupling laser needs to be 
modified from finite to zero to decelerate and stop the input pulse[8, 9]. But when coupling 
laser is very weak, we can no longer treat the coupling laser as classical, in other words, we 
should develop a fully quantized model for EIT. 
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In this paper, both the coupling and probe lasers are treated as quantized. First we give 
a straightforward discussion of the case with no decay rates In section 2.1, and then take the 
decay rates into account in section 2.2. We shall obtain the general form of density matrix 
and susceptibility and analyze several cases in section 2.2: In part(a) we discuss how the fully 
quantized model reduces to the model given by Refs[8, 9] when the mean photon number of 
the coupling laser is large; Part(b) we obtain that the general form of susceptibility and group 
velocity is operators, the expectation value of which can be obtained by act the operators 
on the Fock space and there is fluctuation. We calculate the uncertainty of group velocity 
numerically and give an approximate uncertainty relation between the phase operator of 
coupling laser and the group velocity; Part(c) we discuss the more general case where both 
the probe and coupling lasers are weak and have similar intensities. 



2 The theoretical model 

2.1 The case with no decay rates 

Let us start from the well-known three-level A-type configuration atom (Fig.l) whose energy 
levels assumed to be E a > E c > Ef,. They interact with two quantized fields, probe and 
coupling ones. The two low levels \b> and |c> are coupled to the upper one \a> separately 
and initially the atom is in the ground state 1 6> . The frequency of the coupling laser u 2 = w ac 
and the probe laser oj\ = to a b — Ai, where Ai is the detuning of the probe laser. Here both 
the probe and coupling lasers are quantized. In the interaction picture the Hamiltonian of 
the system is[12]: 

H 1 = Hq + H[ (1) 

Hi = E a \a><a\ + E b \bxb\ + E c \cXc\ 

— h[Ai(\a><a\ + \cXc\)] + uj\a\ai + uj 2 a 2 a 2 (2) 



H[ = HAi(\a><a\ + \c><c\) - h(giai\axb\ + g 2 a 2 \a><c\ + R.C) (3) 

where Oj and a[ (i = 1, 2) are the annihilation and creation operators of the probe (for i = 1) 
and coupling (for i = 2) laser modes respectively, and g\ the coupling constants. Assuming 
the detuning Ai small, we shall find the solution for the dark-state |^o> of the system by 
perturbative approximation. The perturbation with the first order was given in Ref.[ll]: 

i^o> = \4 0) >+\4 1] > 

= —\b,n 1 ,n 2 > ^3— Ai|a,ni - l,n 2 > - — \c,ni - l,n 2 + 1> (4) 

With the energy eigenvalue h-^Ai, i.e., H(\ipo>= h-^Ai\ipQ>. The second-ordered pertur- 
bation can be calculated: 

IVo> = |^> +ivd> +lv>o> 

-A 1 ) \b,ni,n 2 > —r—Ai\a,n 1 -l,n 2 > 



Q 1 40i^ a2 V 

+ -^A? |c,m-l,n 2 + l> (5) 
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where \ni,n 2 > is the usual two-mode Fock basis. Qi = 2g\^fn{ , Q 2 = 2<7 2 \M2 + 1 and 



O = \J$l\ + Q 2 . Suppose the coupling and probe lasers are in a two-mode coherent state 
\a,j3> witha,/? the real for simplicity [12] and the atom is initially in the ground state 
\b>®\a,{3>. If we consider the ideal case in which the decay rates of various levels are 
ignored and initially fii = with Q 2 finite, then proceed to turn £l 2 down while slowly 
turning on. During the course, the state |^o(^)> of the system will evolve adiabatically, 
so from (4) or (5) we can get the density matrix of the system. Further we can calculate the 
susceptibility of the system. If the mean photon number of the coupling laser n a = a 2 and 
the probe laser np = (3 2 are large, i.e., in semiclassical limit, from (5) we have 



P<*W = - 7^T^v? A i + ,o2To2^ A l ( 6 ) 



where 0^ = Qi(n a ,ng) (i = 1,2) are the Rabi frequencies of the coupling and probe lasers, 
with p ab Np ab (ui) = e x(vi)Ei(ui), The susceptibility is given by: 

4N\p ab \ 2 ti 2 lQ\p ah \ 2 N^j 



and 



dwi he (n 2 + n 2 ) 2 he (nj + n 2 f 

The above result is valid for small Ai.The last term of r.h.s of (8) which was not included 
in the Ref.[ll] is negative and it indicates that d ^^ 1 - ) < and then group velocity of the 
probe laser may be greater than the vacuum speed c if we extrapolate (8)to large Ai for the 
abnormal dispersion we meet here. 

The case we discussed above is very ideal. However, in general the decay rates of various 
levels cannot be ignored. In this case, to obtain the susceptibility of the media, we should 
solve the evolution equation of the density matrix. 

2.2 The case with decay rates 

In the Schrodinger picture, the dynamics of the system is described by the interaction Hamil- 
tonian: 

Hi = -N j ^-{hg l aie i ^ {z ' ct) \axb\ 

+hg 2 a 2 e i ^ {z - ct) \a >< c| + H.c) (9) 



where io\ = u ab — Ai, u> 2 = u a c, 9 = Py 2l^F an< ^ ^ ^he quantization volume, N the number 
of atoms in this volume and L its length in z direction. The density matrix of the atom 
system is defined by: 

p(z, t, t ) = PaA z ' *> to)\a><P\ (10) 

a, /3 

where a, (3 = a, b, c, and p a p are the density matrix elements. Make the substitutions: 
Pab = Pabe~ l( - Wab ~ Al ' >t , Pcb = /5 c 6e _i ^ c6_Al )*, and others, p^ u = p^e" 1 ^^ 1 . If the initial state 
of the atom-field system is assumed to be \b> <8| «*,/?>, i.e., &i(t = 0) = 0, fl 2 (t = 0) > 0. 
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In the near-resonance case, very little atoms are populated in the state \a> and the matrix 
elements p~bb and/5 cc varies slowly with t, therefore these two density matrix elements can be, 
respectively, replaced by their initial value p$ and pic which are given by (4) or (5), then 
the evolution equations of the three density matrix elements p a b, p c b and p ca can be written 
in the matrix form: 



R 



-MR + A 



(11) 



where 





' Pab' 




R = 


Pcb 


, M = 




.Pca. 





7i + iAi 

-ig 2 a\e~ lk2 * 



-ig 2 a 2 e ik2Z 
73 + «Ai 
igia\e~ tklZ 





igi&ie 

72 



ik\z 



A 



^P ( bbEi(z) 




(12) 



and E„ 



(m = 1,2), 71, 72 and 73 are the off-diagonal decay rates for 



J m(z) — \f 2 e V 

Pab^ Pca and p c b respectively. Conventionally, both the coupling and probe lasers were treated 
as classical, it should require the coupling laser is much stronger than the probe laser, hence 
only p a b and p c b are needed [4]. However in our case, the coupling and probe laser are both 
quantized, so they may be equally strong, therefore, besides p a b and p c b-, we should consider 
p ca as well. 

When the matrix M is non-singular, the formal solution of equation (11) is given by: 



R = e - Mt R + (1 - e~ mt )M- L A 



Mt\ 



(13) 



where Rq is the initial value of R(t) given by (4) or (5). Let us first make analysis of solutions 
of (13) varying with t: 

i) when t is very small, i.e., e~ Mt ~ 1, so R(t — ► 0) = Rq. However, because of the decay 
rates, only initially R can be given by (4) or (5) and its form will be changed when time 
becomes large. 

ii) when t is large, i.e., e~ Mt ~ 0, then 



R = M' l A 



(14) 



From (14) we see that R(t) will reach a steady value when the time is large enough, i.e., R 
become independent of time. Under such condition we should carefully consider various cases 
as follows: 

(a) If the mean photon number of coupling laser n a = a 2 is large while that of the probe 
laser rip = (5 2 is small, and n a 3> ftp, therefore the influence of the probe laser in equation 

(11) can be ignored and most atoms populated in the ground state, i.e., p$ = l,pic = 0, so 
the matrix M and A are given by 

a ik 2 z 



M 



71 + iAi 

-ig 2 a 2 e~ lk2Z 




-ig 2 a 2 e 
73 + »Ai 







72 





r^£i(z)l 













We then obtain 



Pab 



; Pab 

■ 2h 



( 73 + iAi)£7i(z) 



(71 + iAi)(73 + iAi) + g\a 2 a} 2 



(15) 



(16) 
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Together with p a bNp a b = eoxi^ijE^z), the susceptibility is given by 

ig 2 1 N( l3 + iA 1 ) 



(71 + iAi)(7 3 + iAi) + g%a 2 a\ 



(17) 



which is a formal solution and the susceptibility is an operator . Because n a = a 2 is large 
and the relative fluctuation of the photon number 4^ = ^ is small, we may, as a good 

approximation, replace <7 2 a 2 al, by g\{n a + 1) when calculating the mean value of x( w i) m the 
two-mode coherent state \a,/3>. Noting that f2 2 = g2\/n a + 1, then 



X(^i) 



ui [(71 + iAi)(73 + iAi) + 



(18) 



The above result is familiar [4], which indicates that when the mean photon number of 
the coupling laser is large, our model reduces to the model where the probe laser is quantized 
while the coupling laser is classical [8, 9]. 

(b) If both the probe laser and coupling laser are very weak, but still n a 3> rip, the formal 
solution of the susceptibility has the same form as (17), however, meanwhile we should 
not ignore the relative fluctuation of photon numbers of the coupling laser. The relative 
fluctuation of the susceptibility may also be large, and <7 2 &20 2 can 110 l° n g er be replaced by 
g\ {n a + 1) in the calculation of the mean value of xi^i)- In fact, we have in this case 



X(wi) 



<a, P\x(ui)\a, P> 



E 



^1^(73 + »A a 



a 



2n 2 



U2- 



u! [(71 + *Ai)(ts + iAi) + g\{n 2 + 1)] 



n 2 ! 



(19) 



Noting that x( w i) = Xi^i) + 2X2(^1)1 Xi( u; i) and X2^2) are, respectively, the real and 
imaginary parts of the complex susceptibility and related to the dispersion and absorption: 



Xi^ 1 )=g\Ne- a2 £ 



(73 



A?- 52 2 (n 2 + l))A 1 



a 



2n 2 



n-2 



=0 



and 



X2(^i) 



g 2 Ne~ 



E 



( 7l73 - Af + 5 2 (n 2 + l)) 2 + (71 + 73) 2 A? 



(7i7 3 2 +5 2 (n 2 + l))A 1 



ra 2 ! 



2n 2 



«2 : 



=0 



(7173 - A 2 + 5 2 (n 2 + 1)) + ( 7l + 73 ) 2 A 2 



ra 2 ! 



(20) 



(21) 



Noting that P\ = Axi(k>i)/|xi(wi)| an d P2 = Ax 2 (^i)/|X2(^i)| are the relative fluctuation 
of Xi(^i) and x 2 (u;i) respectively, in Fig.2 and Fig.3, the relative fluctuation of xiO^i) and 
X 2 (^i) are plotted versus the detuning Ai in units of the atomic decay 71 respectively, for 
a 2 = 500, 71 3> 73 and fi 2 = g2\fn a + 1 = 7i/2. It is seen that, around the zero detuning, 
for example, Ai — O.I71. The relative fluctuation of xiO^i) is small (« 4%), while that 
of x 2 (o>i) is large (« 9%). On the other hand, around the detuning Ai = —0.7 the relative 
fluctuation of xi( w i) is large (~ 200%), while that of X2(^i) is small (~ 0.4%). Furthermore, 
the derivative of xii^i) is related to the group velocity for the probe laser pulse through 



V B = C/ 



CLLd 1 



(22) 



from which we can further numerically calculate the accompany fluctuation of the velocity. 
For example, on the zero detuning Ai = 0, we obtain V g = 10.02 m/s, uncertainty AV g = 
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0.45 m/s, and relative fluctuation AV g /V g ~ 4.5%, while on the detuning Ai = O.I671, we 
obtaint^ = 39.45 m/s, uncertainty AV g = 27.90 m/s, and relative fluctuation AV g /V g « 
70.7%. 

The results above shows that the group velocity of the probe laser is not a certainty in 
the fully quantized model. Its uncertainty is a function of detuning Ai. In what follows we 
shall give an approximate uncertainty relation between the phase operator of coupling laser 
and the group velocity, for the two-mode coherent state \a, (3>, as an approximation, we have 

V g « V g (n a , Ai) + d /( !^ 2 Al) ln 2 =n a (n 2 - n«) (23) 

where /(n 2 , Ai) = c/[l + ( 1/2) ^ ( ^7^1^^^ )] , and n 2 = aj«2 is the particle num- 
ber operators of the coupling laser modes and it satisfies the commutation relation[14]: 

[n 2 ,cos0] = — zsin</> (24) 

where cos <j> is the phase operator of the coupling laser, then 

[y 9 ,cos0] = -iF(n Q ,Ai)sin<£ (25) 

where F(n a , A\) = 9 ^q^^ \ n2 =n a , from the uncertainty principle, we have 



(AV^)(Acos0) > ^(na.AiXsin^l (20) 



which turns out that the uncertainty of the group velocity is the function of n a and Ai. 

As we have known that to decelerate and stop the input pulse, the strength of the coupling 
laser need to be modified from finite to zero [8, 9, 13], but when fi 2 is very small, we should 
treat the coupling laser as quantized. It is noticeable that if initially the coupling laser is 
much strong than the probe one, the coupling laser will be much stronger than the probe 
laser at all times(see Ref.[8]), which satisfies the condition of part (b) we discussed above. 

(c) When both the probe and coupling lasers are very weak and have similar intensities, 

i.e., n a ~ np, from (4) we have p$ = ^§ and p^c = hence from the equation (11) the 
form of the density matrix element p a b can be followed: 

il&EkWPbb [(73 + iAx) 72 + gfa^l 



2% 

Pab 



(71 + iAi)(73 + iAi)72 + 72 3 2 a 2 4 + ( 7 i + iA^gjaial 
and 



(27) 



x(.^i) = — r :— 71 7+ — ; — - — . ; „. .+i ( 28 ) 



(71 + «Ai)(73 + iAi) 72 + 7 2 c/|a 2 4 + (71 + i&i)glaia\ 



From the above result we find that the susceptibility depends on the coupling laser as well as 
the probe laser when both of them are equally strong, which is similar to the result of KCW 
in Ref[ll], but here the susceptibility is an operator and its value can be calculated on the 
Fock space just as the analysis in part(b). 



6 



3 Conclusion 



We have developed a fully quantized model for EIT in which the decay rates are taken into 
account. In this model, to solve the evolution equation for density matrix, we separate the 
atom-system from the photon-system described by Fock-states and only calculate the density 
matrix of atom-system. By this means, the general form of susceptibility and group velocity 
of the probe laser we obtained are operators concerned with particle number operators of the 
probe and coupling laser modes. Their expectation value and fluctuation can be calculated 
on the Fock space. We have calculated the uncertainty of the group velocity numerically and 
give an approximate uncertainty relation between the phase operator of the coupling laser 
and the group velocity. When both the probe and coupling lasers are weak and have similar 
intensities, we find the susceptibility depends on the coupling laser as well as the probe laser. 
Considering the decay rates of various levels, we can make analysis of absorption of probe 
laser near resonance and calculate the fluctuation in both absorption and dispersion. We 
also discuss how the fully quantized model reduces to a semiclassical model when the mean 
photon numbers of the coupling laser is large. 

This work is in part supported by NSF of China. 
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Figure 3: the relative fluctuation of X2(^i), for a 2 = 500, 71 S> 73 and ^2 = 92^/n a + 1 = 
7i/2. 
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